A particle-number-conserving solution to the generalized pairing 

problem 



Feng Pan*, J. P. Draayer, and W. E. Ormand 

Department of Physics and Astronomy, 202 Nicholson Hall, 
Louisiana State University, Baton Rouge, LA 70803-4001 



Abstract 

An exact, number-conserving solution to the generalized, orbit-dependent 
pairing problem is derived by introducing an infinite-dimensional algebra. 
A method for obtaining eigenvalues and eigenvectors of the corresponding 
Hamiltonian is also given. The relevance of the orbit-dependent pairing solu- 
tion is demonstrated by comparing predictions of the model with shell-model 
calculations. 
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Pairing has long been considered an important interaction in physics. The concept can 
be traced back to the seniority scheme introduced by Racah in atomic physics ||T|. Its 
physical significance was first realized in the study of superconductivity |2j . Following the 
suggestions of Bohr, Mottelson, and Pines ||, the first detailed application of pairing in 
nuclei was made by Belyaev [|J. The concept has since been applied to other phenomena: 
high T c superconductivity [HH, applications using the Hubbard model 0, and pairing 
phenomena in liquids f| and metal clusters f|. 

BCS methods have yielded major successes in studies of superconductivity. When ap- 
plied to nuclei, however, one must remember that the number of spectroscopically active 
nucleons is typically too small {n ~ 10) to support underlying assumptions of the theory, 
namely Sn/n is not negligible, and as a consequence particle-number nonconservation effects 
can enter and give rise to spurious states, nonorthogonal solutions, etc. Although this chal- 
lenge can be partially remedied by making use of particle-number projection techniques [[TO 



the elegance and simplicity of the BCS method are then compromised. Secondly, an essential 
feature of pairing correlations are even-odd differences, which are driven mainly by Pauli 
blocking. It is difficult to treat these differences in the BCS formalism because different 
quasi-particle bases must be introduced for different blocked levels. After an investigation 
into the accuracy of the BCS approximation, Kerman and Lawson suggested that an exact 
diagonalization of the pairing Hamiltonian is necessary Jl0| . Based on these observations, a 



particle-number-conserving method for treating the pairing problem in well-deformed nuclei 



was put forward (TT| . The method uses a configuration-energy truncation scheme and takes 
the strength of the pairing interaction to be the same for all orbitals. Unfortunately, because 
of the deformation, each orbital can only accommodate a single pair of particles and this 
limits the applicability of the theory. Complementary mean-field techniques with approxi- 
mate Lipkin-Nogami number projection and non-trivial pairing interactions have also 



been found to be a useful means for handling pairing correlations in such systems fll^HlB 
However, exact shell-model solutions, when feasible, even if only approximate, are the best 
way to probe the true nature of many-body correlation effects. 



The purpose of this contribution is to give an exact, particle-number-conserving solution 
to the generalized, orbit-dependent pairing problem. The method can be applied to any 
number of pairs in any model space, even to cases where exact diagonalization is computa- 
tionally prohibitive. 

The generalized pairing Hamiltonian for spherical nuclei can be written as 

H = e 3 a ]m a jm ~ \G\ CjfS + {j)S~{f), (1) 

jm jj' 

where the 6j are single-particle energies and S ± (j) and S°(j) are the pairing operators for a 
single- j shell defined by 

m>0 

^ U) = ( — y ma j-m a jm, 

m>0 

1 1 - 

S °U) = o £ ( a ] m ajm + a]_ m aj_ m - 1) = -(Nj - Qj), (2) 

where Q 3 ■ = j + 1/2 is the maximum number of pairs in the j-th shell, Nj is the j-th shell 
particle number operator, and the cjji measure the orbit-orbit pairing strength. Note that 
Hamiltonians with degenerate single-particle energies (ej = e for all j-orbitals) and separable 
pairing strengths (cjj> = c*Cj>) are special cases of the general theory. 

In general, for N pairs, Eq. ([[]) can be diagonalized with bases states that are products 
of the single- j shell pairing wave functions: 

\N) = Eflfefc-^ k ^S? 2 k *S+ k * ..-St fe *|0), (3) 

where the summation is restricted so that 

Eh = N, (4) 

i=l 

and the i^fofa— fc are expansion coefficients to be determined. Here, |0) is the pairing vacuum 
state that satisfying the condition 

S~\0) = for all j. (5) 
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The dimensionality of the Hamiltonian matrix in this basis increases rapidly with increasing 
N and the number of shells p. Due to the Pauli Principle, it is less than or equal to the 
dimension of the irreducible representation (irrep) [NO] of the unitary group U(p), 

dim < (p+^-y. ( 6 ) 

The equal sign in Eq. (Jg) holds when all the single- j shell pairing wave functions in the 
summation of Eq. (||]) are Pauli allowed. From Eq. (g), it is clear that the problem quickly 
becomes intractable because there are no analytical expressions or recursion relations for 
determining the Bk^-kp coefficients. 



As in the quasi-spin case |p.7|| , we consider a slightly simpler Hamiltonian 

H = a jm a jm - | G | So" So > (?) 

jm 

that is, one with degenerate single-particle energies and a separable pairing interaction, 

# = £#^0"), s =J2c 3 s-(j) (8) 

3 j 

with the coefficients q satisfying the condition 

£M 2 = i- (9) 



This defines generalized pairing as proposed by Talmi |]18| . The separability assumption, 
though strong, is physically motivated as it links the pair-pair interaction strength to the 
individual pair formation probability. In the notation of Eq. (|]), Cjy = c*cy with \cj\ 2 
giving the percentage of single-j shell pairing in the Hamiltonian. In what follows, the Cj 
are taken to be real. Because the total number of particles is a good quantum number in 
a number- conserving theory, the single-particle term in Eq. (0) is a constant and can be 
dropped without loss. Hence, Eq. (|7|) reduces to 

H = -\G\S+S . (10) 

To diagonalize this Hamiltonian, consider an algebra generated by 
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S± = Ec* m+1 S ± (j)- (11) 
i 

It is easy to show that these generators satisfy the following commutation relations: 

[S m , S n ] = 2S m+n+1 , 

[S m , S n ] = ±S m+n . (12) 

Therefore, the {S^, \l — 0, +, — ; m — 0, ±1, ±2, • ■ •} form an infinite-dimensional algebra, 
one that differs only slightly from a general Lie-algebra of the affine type without central 
extension. 

The unique lowest-weight state of this algebra is simply the product of the single- j shell 
pairing vacua with arbitrary seniority quantum numbers. Therefore, it suffices to consider 
the total seniority zero case. The lowest-weight state satisfies 

S m \0} = 0; m = 0, ±1, ±2, (13) 

and 

^|0) = -iEl^f m ^|0)=A m |0). (14) 
z i 

Furthermore, it can be proven that the eigenvectors of H for any iV and non-zero energy 
eigenvalue can be written as 

\N)=AfS+S+S+ 2 ...S+ N JO), (15) 

where M is a normalization constant and 

^=E^ (J ). (16) 



To obtain the variables {xf, i = 1, 2, • • • , iV — 1}, Eq. (15) can be expanded in terms of 



Xi around Xi = 0, 



\N) = ATEzFx? ■ ■■x^S+S^ ■ ■ -5+ _J0), (17) 



where the are the Fourier-Laurent coefficients in the expansion of S£ , namely 

a t = -Li' b "? s i- (18) 

Using Eq. ( ]T7| ) and the commutation relations of Eq. flI2|), it is easily shown that the Xj, 
with i — 1, 2, • • • , N — 1, satisfy the relations 

-;E" A 1 2 =- + E^— . < = 1, 2,.-., JV-1, (19) 

with 

Ss =1> (20) 

where 

W = x,/a, a = -^-^-, /» = J5/(-|G|). (21) 

We note that although these relations were derived for Xi ~ 0, they are nonetheless valid in 



the entire complex plane except at the singularities in Eqs. (19) and (|20|) . Therefore, the co- 
efficients Xi (i=l, 2,- ■ -, N-l) and eigenvalues of the pairing energy E ^ are simultaneously 
determined by the system of equations Eqs. (|l!J) and (|20|). 
Similarly, the eigenvectors for E = can be expressed as 

\N, 0>= ^5+52,-.. 5^10). (22) 

Using the same technique as in the E ^ case, it can be shown that in this case the Xi with 
i — 1, 2, • • • , N, are determined by the following set of equations 



c 2 



for = 1, (23) 



|E%— rn = Er4^-> ^ = 1, 2, iV, for iV>2. (24) 

At this point, we reiterate that our solution is for a separable, orbit-dependent pairing 
interaction as defined by Eq. (|7|) , and not for the generalized pairing Hamiltonian of Eq. ([[]) . 
If the separability assumption is valid, the solutions should be good approximations to the 
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more general theory. The power of the method lies in the fact that it gives the exact, 
particle-number-conserving solution of the pairing problem for any number of pairs in any 
model space, even for cases where full matrix diagonalizations cannot be carried out. Next 
we show, by example, that the separability assumption is consistent with "realistic" shell- 
model results and hence that theory can be used with confidence to study many-particle 
correlation phenomena in pair-rich systems. 

As a test of the method, comparisons with shell-model calculations for like particles 
interacting through a "realistic" generalized pairing interaction were carried out. The shell- 
model calculations were performed within the configuration space defined by the nuclear 
<is-shell (0d 5 / 2 , 0^3/2, and lsi/2 orbitals). The generalized pairing Hamiltonian of Eq. ([[]) 
was defined using the J = two-body matrix elements of the universal cfe-shell Hamiltonian 
of Wildenthal ||19|| . Since the number of J = states for any two-particle system is small 



(three for the ds-shell), it is reasonable to define the generalized pair in Eq. (H) as the wave 
function of the shell-model ground state. For our ds-shell example this means the coefficients 
are: q 5/2 = 0.570243, q 3/2 = 0.626089, and c Sl/2 = 0.531823. The strength of the interaction 
was then adjusted to reproduce the two-particle, ground-state energy, that is, \G\ = 3.12778. 
With the Hamiltonian so defined, solutions were obtained for N = 1, 2, and 3 pairs of like 
particles ||20|| . Shown in Fig. |l| is a comparison between the eigenenergies obtained from the 
shell model (SM) and the separable, generalized pairing method. For reference, the quasi- 
spin (QSA) approximation gives -6.389, -10.684, and -12.788, respectively, for the N = 1, 



2, and 3 pair cases |2T|. A further measure of how well the generalized pairing solutions 
correspond to the shell-model results is illustrated in Table | where the percent overlap 
between the shell-model and generalized pairing wave functions are given. In all cases, the 
wave functions obtained with the method presented here have better than a 94% overlap 
(96% for the ground state) with the shell-model states, which indicates that the assumption 
of separability is sound. 

An exact, particle-number-conserving solution to the generalized, orbit-dependent pair- 
ing problem has been derived. A key feature of the derivation is the use of an infinite- 



dimensional algebra, which is new in nuclear physics. Because the theory is number- 
conserving, it goes well-beyond what has been done even for the case of equal pairing 
strengths which includes to quasi-spin as a special limit. The relevance of the theory was 
shown by comparing eigenvalues and eigenvectors with "realistic" shell-model results. 
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TABLES 

TABLE I. Percent overlap between the eigenstates obtained with the separable, generalized 
pairing method and the shell model as a function of the number of pairs, N. 



eigenstate \ N 1 2 3 

1 96.2 99.4 97.8 

2 94.7 97.6 97.0 

3 99.4 97.2 94.3 

4 98.8 95.8 

5 99.2 96.8 

6 97.0 
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FIGURES 

FIG. 1. Comparison between the spectra obtained from the shell model (SM) and the separable, 
generalized pairing (SGP) method as a function of the number of pairs, N. The dashed lines indicate 
the correspondance between the levels in the two models. 
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